.mage . -pp 



'{ 

for(i-0; i<m_nuraPixels ; 

{ . 

if (show_progress && i%100==0) load_progress . SetPos ( (99*i > /m_numPixels) ; 
p=GetPixel (i) ; if (p>=color_map_size) p=*color_map_size-l ; 
SetPixel (i, color_map [p] ) ; 

} 

} 

return true; 

} 

} 

bool Image :: SetPalette (int px, int py, int pxmax, int pymax) // set palette with respect to the mouse posi 
ion 

{ . 

int of f set= ( ( { long) nwnaxPixe lvalue/ 3 ) * ( 2*px-pxmax) ) /pxmax; 
double x= (2 . 0*py) /pymax; 
. x=l+2* (x-1) * * (x-1) ; . 
return m_pPal->SetPalette (of f set , x) ; 

} 



* Plot DIB pixels 

bool -Image :: DisplayDIB (CDC v pDC, CRect crectScreenArea, CRect erect ImageArea, 

CPoint- pScroll, bool optimize) 

{ 

U INT palUsage; 

if (m_DisplayFailed) return false; 

if {m_pPal->p_active m_pPal->LoadPalecte (pDC, m_RGB) ) 

{ • 

palUsage=DI3_PAL_CGL0RS; 

} 

else palUsage=DIB_RGB_COLORS; 

/* Optimize drawing regions */ 
if (optimize) 
{ 

CRect rcClip, rcDraw, rcDIB; 

pDC->GetClipBox (rcClip) ; 
rcClip. MormalizeRect () ; 

if (rcClip. IsRectEmpty () ) return true; 
rcClip. InflateRect (2, 2) ; 

rcDraw. intersectRect (rcClip, crectScreenArea) ; 
rcDraw. MormalizeRect () ; 

if (rcDraw. IsRectEmpty {) ) return true; 

rcDIB=m_smS.Screen_to_Image (rcDraw) ; 
if (rcDIB. IsRectEmpty () ) return true; 

crectScreenArea . CopyRect \ rcDraw) ; 

if (theApp. app_Metheus) crectScreenArea . Of fsetRect (-pScroll) ; 
erect ImageArea .CopyRect ( rcDIB}'; 

} ■ - 

/+ Correct image area */ 

if (crectlmageArea. left<0) erect ImageArea .Of fsetRect ( -crectlmageArea . lef t, 0) ; 
else if (crectlmageArea. r.ight>ra_Width) 

crectlmageArea. Of fsetRect (m_Width-crectImageArea. right, 0).; 
if (crectlmageArea. top <0> crectlmageArea .Of fsetRect (0, -erect ImageArea. top) ; 
else if (crectlmageArea . bottom>ra_Height ) 

crectlmageArea .Of fsetRect (0, m_Height-crectImageArea . bottom* ;. 

/* Set screen (destination) and bitmap (source) areas */ 

long xDest^crectScreehArea.TopLeft (> .x; long yDest=crectScreenArea . TopLef t (). y; 
• long wDest^crectScreenArea. Width { ) ; long hDest=crectScreenArea . Height () ; 

long xBmp=crectImageArea. TopLef t ( ) .x; 
long yBmp; 

if (theApp. app_Metheus) yBmp=crect ImageArea . TopLef t (). y; ■ EXHIBIT 

f se h A 

yBmp=this->GetHeight () -crectlmageArea. TopLeft () . y-crect ImageArea. Height () ; J| 
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long wBmp^crectlmageArea.Widti, . t ; long hBmp=crect ImageArea . Heig. 

switch (ra_Flip) 
{ 

case 1: // horizontal 

xDest += wDest-1; xBmp=this->GetWidth ( ) -xBmp-wBmp; 
wDest = -wDest; break; 
case 2: // vertical 

yDest += hDest-1; yBmp=this->GetHeight ( ) -yBmp-hBmp; 
hDest =» -hDest; breaks- 
case 3: // 130 

xDest wDest-1; xBmp=this ->Get Width () -xBmp-wBmp; 
wDest « -wDest; yDest +=• hDest-1; 

yBmp=thi s->Get Height () -yBmp-hBmp; hDest = -hDest; break; 

} 

/* Display as DIB */ 

if (wDest==0 II hDest==0 II wBmp ==0 II h5mp==0) return true; 
if ( ! theApD . app_Metheus } 
{ 

m_DisplayFailed= {StretchDIBits ( pDC->GetSaf eHdc ( ) , 

xDest, yDest, wDest, hDest, xBmp, yBmp, wBmp, hBmp, 
m_Pixels, m_bmpInfo, palUsage, SRCCOPY) <=0) ; 

} 

else // Metheus 
I 

m_DisplayFailed=true; 
if (m_apdate3itmap) 
{ 

if (m_Bitmap) 

r-- m , . MetheusDeleteCompatibleBicmap ( pDC->Get5af eHdc f ) , m_Bitmap) ; 

m_Bitmap = MetheusCrea teCompatibleBitmap t pDC->GetSaf eHdc t ) , GetWidth ( ) , GetHeight ( ) i , 
if {MetheusLoadlmageFromData ipDC->GetSafeHdc ( ) , m_3itmap, 

theApp. app_DynamicPalet teStart , m_Pixels, GetWidth 0/ GetHeight (), 
GetWidth ( ) *m_3ytes_per_ Pixel, 8*m_3ytes_per_Pixel, 
0, 0, GetWidth (i, GetHeight i), 0,0)«FALSE) 

{ 

AfxMessageBox ( "Cannot create image bitmap", MB_0K I MB_ICONEXCLAMATION) ; 
return false; 

} 

m_Di splay Failed= 

(MetheusStretchBltlmage ;pDC->Get5af eHdc ( } , NULL, 

xDest, yDest, wDest, hDest, 
m_Bitmap, 

xBmp, yBmp, wBmp, hBmp, SRCCOPY) == FALSE) ; 
#* 

} 

if (m_DisplayFailed) // Display error 
{ 

DWORD ecode=GetLastEr < ror ( ) ; 
CString estr; ' " 

estr . Format ( " GDI error code: *ld\n Please reload the image", ecode) ; 
AfxMessageBox (estr, MB_OK I MB_ICONEXCLAMATION) ; 

} 

m_Upda t eB i tmap= false; 
return { !m_DisplayFa^lg4i*»fc.- 

bool Image: :DisplayDIB {CDC *pDC) 

return DisplayDIB (pDC,m_smS.crScreen,m_smS.cr Image, CPcint(0,0), true); 

bool Image: : DisplayDIB (CDC ^pDC, CPoint pScroil) 

return DisplayDIB (pDC, m_smS . crScreen, m_smS . crlmage, pScroil, true) ; 



■# 

* Get image width 

+ + + * + + + + + + + + * + * + + + *^ + ** + + + * + + + * + + ** + * + **+ + + + + »^ 

int Image: : GetWidth () 

{ 

return m_Width; 

} 
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si +- t*l. 09861229; s2 +- t; 
// d-sqrt(lO) 



p=Get Luminance (x-1, y-3) 
p=Get Luminance (x-1, y+3) 
p=Get Luminance (x+1, y-3) 
p^GetLuminance (x+1, y+3 J 
p=^GetLuminance (x-3, y-1) 
p=Get Luminance (x-3, y+1) 
p=GetLuminance {x+3, y-1) 
p=GetLuminance (x+3, y+1) 



pmax=pmin=p ; 



i f ( p>pmax ) pmax=p ; 

if(p>pmax) pmax=p; 

i f ( p>pmax ) pmax^p ; 

if (p>pmax) pmax-p; 

if(p>pmax) pmax=p; 

if(p>pmax) pmax=p; 

i f ( p>pmax ) pmax=p ; 



else if(p<pmin) prain=p; 

else if(p<pmin) pmin«p; 

else if(p<pmin) pmin=p; 

else if (p<pmin) pmin=p; 

else if (p<pmin) pmin=p; 

else if (p<pmin) pmin=p; 

else if(p<pmin) pmin^p; 



t=log (max (pmax-pmin, 1 ) ) ; 
si +- t*l. 15129255; s2 += t; 
Sl *= 0.14235714; s2 » s2*0 . 10477684; 

//Find fractal dimension as slope of the Hurst line 
//t=max(6. 386491206* (sl-s2) , 0.0) ; // fractal dimension 
t^max(600*m_TR_scaie + (sl-s2) , 0. 0) ; // fractal dimension 
return ( (long) (t ) ) ; 



Apply a pixel neighborhood function (with code mask_type) to the image 



bool Image : :TR_PixelNeighboorhood (char mask_type, bool show_progress) 
{ 

int rad, i, j, jl, top, bot; 
long p, pmax, pcovint; 
double pavg; 

// Set pointer to the masking function 
long ( Image :: *maskF) (const int x, const int y) = 0; 
switch (mask_tvpe) 
{ 



case 


*g' : 


maskF= 


=TR 


_Gauss ; 


rad=l; 


break; 


// 


gaussian denoising 


case 


' a 1 : 


maskF= 




Smooth; 


rad=l; 


break; 


// 


average smoothing 


case 


•s' : 


maskF= 


=tr" 


"Sharp; 


rad=l; 


break; 


// 


sharpening 


case 


'e' : 


maskF= 


=tr~ 


"Sobel; 


rad=l; 


break; 


// 


edge detector 


case 


' f * : 


maskF= 


*tr~ 


Fractal; 


rad=3; 


break; 


// 


edge detector 



default: return false; // invalid mask type 
} 

// Create temporary buffer 

int w=this->GetWidth ( ) ; 

int h=this->GetHeight () ; 

long * { *buf ) =new long* [ rad+1 ] ; 

if { !buf ) 

{ 

Af xMessageBox ( "Low memory, cannot transform"); 
return false; 



for(i=0; i<=rad; i++) 
{ 

buf[i] = new long[w]; 

if (buf [i]==0) 

{ 

Af xMessageBox ( "Low memory, cannot transform"); 
fos(j=0; j<i; ( if(buf[jj) delete [J buf[j]; } 

delete [] buf; 
return false; 
} // out of memory 



// Display progress, control in the main frame status bar 
CProgressCtrl load_progressr 
if (show_progress) 
{ 

CreateProgressControl (load_progress, "Transforming ..."); 
load_progress . SetPos ( 5) ; 



// Estimate max transformed pixel value 
m_TR_scale=>l . 0; 

p=pmax= (this->*maskF) (rad+1, rad+1) ; 

pcount=0; ■ pavg=0.0; 

int dw = w>">4; if (dw<l) dw»l; 

int dh = h»4; if (dh<l) dh=l; 

for(i=rad+l; i<w-rad; i += dw) 

{ 



+ + + + # + + + + + + + 




★ 
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for(j^rad+l; j<h-rad; j ih) 
{ 

p=(this->*maskF) (i, j) ; 
pcount ++; 
pavg +» p; 
if(p>praax) pmax=p; 

) 

) 

if(pmax>0 && pcount>l) 
{ 

pavg /= pcount; 

pmax=min (pmax, (long) (2* pavg) ) ; 

m_TR_sca le^O. 6* (double) (m_maxPixe lvalue ) /pmax; 

} 

if (show_progress) load_progress . SetPos (10) ; 

// Apply (2*rad-l) * (2*rad-l) masking operator 

bot=-l; 

top=rad-l ; 

for(j=rad; j<~h; 

{ 

if (show_progress j%50==0) load_progress . SetPos ( 10+ (90* j ) /h) ; 

jl=j-rad-l; 

if{jl>-rad) 

{ 

for(i=rad; i<w-rad; i++) SetPixel (i, j L, buf [bot ] [i] ) ; 

} 

bot = (bot+1) % (rad+1) ; 
top' = (top+1) t (rad+1) ; 
if(j<h-rad) 
{ 

for{i=rad; i<w-rad;' i++) buf [ top] f i ] ~ < this->*maskF) (i , j ) ; 

} 

} 

// Clean up 

.for(j=0; jorad; { if (buf [j]) delete [] buf[j]; } 

delete [] buf; 
Beep(500, 50) ; 
return true; 



^ + + + + + + * + + + * + ^ + + + + + + + + + * + + + + * + ^ 

* Inverts the image 

* + + * + + +* + + + + + + + ++ ++ + *+ + + * + + + + + + + * + + + ±+ * + + + + * + + + ** + ±*-k-*irir + + *'+ + + + + + + + + + / 

bool Image : : TR_Negate ( ) ^ 
{ 

if (m_pPal->p_active) m_pPal->Negate ( } ; 
else 
■ { 

// Display progress control in the main frame status bar 
CProgressCtrl load_progress ; 

CreateProgressControl (load_progress, "Changing to negative image ..."); 
// Invert pixels values 
if (m_RGB) 

{ ' 
BYTE r, g, be- 
long p; ... 
fordong i=0; i<m_numPixels; i++) • 

( ' 
if(i%1000~0) load_progress. SetPos ( (99*i> /m_numPixels) ; 
p=GetPixel (i) ; 

r=m_max Pixel Value -Get Rvalue (p) ; 
g=m_maxPixelValue-GetGValue (p) ; 
b=ra_maxPixelValue-GetSValue (p) ; 
SetPixel (i, RGB (r,g,b) ) ; 

} 

} 

else 
( 
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for (long i~0; i<m_nu: els; 
{ 

if (i%1000 M 0) load_progress .SetPos ( {99*i j /nwiumPixels) ; 
SetPixel (i,m_maxPixelValue-Get Pixel (i) ) ; 

} 

) 

} 

Beep{500, 100) ; 
return true; 



* Copy pixels from current to safe buffer (on true) . 

* or vice versa (on false) 

void Image :: ResetPixeis (bool current_to_saf e) 
{ 

if (current_to_safe) // current -> safe 
{ 

memcpy (m_Saf e_Pixels , m_Pixels , m_numPixel Bytes ) ; 

} 

else // safe -> current 
{ 

memcpy (m_Pixeis , m_Saf e_Pixels , m_numPixelBytes ) ; 
m_UpdateBitmap=true; 



+ 

* Histogramm stretch from [amin,amax} to [bmin, bmax'] color range. 

* R,G and B components are stretched together 



bool Image : :TR_HistStretch (int amin, int amax, int bmin, int bmax, 

bodl show_progress) 

{ 

■ long i, cmax, p; 

. /* Create color map +/ 
if (m_pPal->p_active) cmax=m_pPal->p_Size; 
else Get_Pixel_minmax fi, cmax) ; 
cmax++; 

long* color_map=new long [cmax]; 
• if ( !color_map) 
{ 

AfxMessageBox ( "Low memory, cannot perform this transform"); 
return false; 
} // out of memory 

■ /* Set color map parameters +/ 
if (amin<0) amin=0; 
if{bmin<0) bmin=0;- 

if {amax>m_maxPixelValue) amax^n^maxPixe lvalue ; 

if (bmax>m_max Pixel Value) bmax^n^maxPixelValue; 

if (amin>=amax I I. bmin>=bmax) return false; // invalid. map 

if (amin==bmin && amax M bmax): return true; // no stretch needed 

/* Fill the color map^ V 
long da=amax-amin; 
long db=bmax-bmin; 
for(i=0; Kcmax; i++) 
'{ 

p=bmin+(db* (i-amin) ) /da; 

if (p<bmin) p=bmin; else if (p>bmaxi p=bmax; 
color_map[i] =p; • 

> .* 

SetPalette (color_map, cmax, show_progress) ; 

delete [] color_map; 
return true; 

} 
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* Histograms stretch from (percent) % median neighborhood 

* to .the maximal [0, m_raaxPixe lvalue] range 
•* 

#+++ * ++ *+ + * #+ * +#+++ * + ****++#++**+++*++++++***+*^ * *'/ 

bool Image: :TR_HistStretch (BYTE percent, bool show_progress) 

{ 

long i, amin, amax, amed, p; 

/* Validation */ 

if (percent>=100) percent=99; 

if (percent<0) return" true; 

/* Find pixel statistics */ 
Get_Pixel_minmax tamin, amax) ; 

/* Initialise image histogram */ 

long cmax=amax+l; 

int* hist=new int[cmaxj; 

if(!hist) 

{ 

Af xMessageBox ( "Low memory, cannot perform this transform"); 

return false; 
} // out of memory 
for(i=0; i<cmax; i++) hist [i 3=0; 

/* Estimate image histogram */ ' 
int hmax=200Q0; 

int di= max (l,m_numPixels/10000) ; 

if(m_RGB) // Color image 

for(i=0; i<m_numPixel3vtes; i di; 
{ 

p=m_Pixels [i] ; 

if (hist [p] <hmax) hist[p] di; 

} 

} 

else // Grevscaie image 

{ ' 
for(i=0; i<m_numPixels ; i += di) 

{ 

p=GetPixel ( i) ; 

if (hist [p] <hmax) hist[p] += di; 

} 

} 

if (percent==0) return TR_HistStretch (amiri, amax, 0, 

mjnaxPixelValue, show_progress) ; // simple stretch 

/+ Find histogram color average V 
double ptot=0, tot=0; 
for(i=0; Kcmax; i++) 
{ 

ptot +- ( (double) i) *hist [i] ; 
tot += hist [i ] ; 

> 

amed = (long) (0 . 5+ptot/tot) ; 

/* Find new intensity range to preserve */ 

double keep_max= (100-percent) *tot/100; // number of pixels to keep 

double keep-hist [amed] ; 

long bmin=amed; long bmax=amed; ' 

do // do at least once to* guarantee bmin<bmax 

{ . ' 
if (bmin>amin) 

( 

bmin--; 
■ keep hist [bmin]; 

} 

if (bmax<amax) 
£ 

bmax++; 

keep += hist[bmax]; 

} 

} while {keep<=keep_max) ; 
delete [] hist; 

if ( { (bmin!=amin) I I (bmax!=amax) ) (bmin<bmax) ) 
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return TR_HistStretch <bi bxoax, 0, m_maxPixel Value, show ^progress) 
else return false; 

} 



f + + + + + + +* + + + + + +* + + + + + + + + 4+ + + + + + + + + + + + + + + + + + + + + + + + * + + + * + + *+ + ++ *4 + + + + + + + * + 1r + + + + +1r + + -+ + 4 
+ 

* Histogramm equalization 
+ 

+ + + + * + + + + + + + + + ** + + #****** + ** + + * + ■*■"* + ** + ** + ■•■*** + *♦.#•*# + *** 

bool Image: :TR_HistEqualize (bool 3how_progress) 

{ 

long i, amin, amax, p; 

/* Find pixel statistics V 
Get_Pixel_minmax (amin, amax) ; • 
Beep{300, 100) ; 

/* Initialize image histogram, also used as color map */ 

long cmax=amax+l; 

long* hist=new long[cmax]; 

if ( !hist) 

( 

AfxMessageBox ( "Low memory, cannot perform this transform"); 

return false; 
} // out of memory 
for(i=0; Kcmax; i++) hist[i]=0; 

/* Compute image histogram */ 
long hma;<=2 000000; 

int di= ma>: ( 1 , nwiumPixels/ 10000) ; 

if(m_RGB) // Color image 
{ 

for(i=0; i<m_numFixelBytes ; i di) 
{ 

p=m_?ixels [ i ] ; 

if (hist [p] <hmax) hist[p] += di; 

} 

} 

else // Greyscale image 
{ 

for(i=0; i<m_numPixels ; i += di) 
( 

p=GetPixel (i) ; 

if (hist [p]<hmax) hist(p] += di; 

} 

} 

/* Integrate the histogram '/ 
double htotal=0; 

for(i=0; Kcmax; i++) htotal += hist[i]; 

if(htotaKl) return false; // did we have negative pixels or empty image ? 
if (hist [0]<htotal-l) { htotal -= hist[0]; hist[0]=0; } 

/* Fill the color map */ 
double hcum=0; 
for(i=0; Kcmax; i++) 
{ 

hcum += hist[i]; // update cumulative hist 
hist [i] = (long) ( (m_maxPixelValue*hcum) /htotal) ; 

} 

Remap the pixel data */ 
SetPalette (hist, cmax, show_progress) ; . 

/* Clean up */ 
delete [J hist; 
return true; 



* 

* Image Gamma correction 
+ 

bool Image : : TR_GammaCorrect ion (double gamma) 
( 
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#if ! defined (AFX_LUPA_H 24CE8B5 08 F 1 1D2_958F_000000000000 INCLUDED^ 

#define AFX_LUPA_H 24CE8B94_7D8F_11D2~958F_000000000000 INCLUDED^ 

#if _MSC_VER >= 1000 
#pragma once 

#endif // _MSC_VER >= 1000 
// Lupa.h : header file 
// ' 

# include n Image. h" 

///////////////////////////////////////////////////////////////////////////// 
// Lupa dialog 

class Lupa : public CDialog 
{ 

// Construction 
public: 

LupatCWnd* pParent * NULL); 
-Lupa ( ) ; 

// Dialog Data 

/ / { { AFX_DATA ( Lupa ) 

en urn { IDD ~ I DD_MAGN I FY_DI ALOG }; 
private : 

BOOL l_optimize; 
long l_height; 
long l_width; 
public: 

double l_zoom; 
//} }AFX_DATA 



// Overrides 

// ClassWizard generated virtual function overrides 

/ / { { AFX_VI RTUAL ( Lupa } 

protected: 

virtual void DoDataExchange (CDataExchange* pDX) ; // DDX/DDV support 
//} }AFX_ VI RTUAL 

// Implementation 
public: 

bool l_active; 

CSize l_scnSize; 

void Initialize (CSize csScn, double scrn_zoom, double lupa_zoom) ; 

void Reset_l_DC (CDC *pDC, CRect& scrolled_client) ; 

void MoveiCPoints a, CPoint& rel, Image* pBmp, CDC* pDC) ; 

bool Resize (CDC* pDC, CPoint center, CPoint vertex); 

bool Resize (CDC • *pDC) ; 

CString toStringO; 

inline CRect Lupa :: SetlmgRect (CPoint & cpl) 
{ 

CRect r( CPoint (cpl .x- (l_imgSize. cx»l ) , cpl . y- ( l_iragSize . cy»l ) ) , l_imgSize ); 
return r; 

}; 

inline CRect 3etScrnRect (CPoint & cpS) 
{ 

• return CRect (CPoint ( cpS . x- ( this->l_scnSize . cx»l J , cpS . y- ( this->l_scnSize . cy»l ) ) , 
this->l_scnSize) ; 

}; 

protected: 

// Generated message map functions 
//{ (AFX_MSG(Lupa)- 
virtual BOOL OnlnitDialog ( ) ; 
afx_msg void OnChangeMagnif yWidthEdit () ; 
afx_msg void OnChangeMagnif yHeightEdit () ; 
afx_msg void OnChangeMagnif yZoomEdit () ; 
afx_msg void OnChangeMagnif yOptimize () ; 
//} }AFX_MSG 
DECLARE _MESSAGE_MAP ( ) 
private: 

bool l_preactive; 

double l_img_zoom; 

CSize l_imgSize; 

CRect l_scnRect, l_dragRect; 

HBITMAP 1_DIB; 

BYTE* l_Data; 

CDC* 1_DC; 

void Draw(CDC *pDC) ; 
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bool Update2 {CRectS a, CRe-.<ib); 

}; 

// { {AFX_INSERT_LOCATION} } 

// Microsoft Developer Studio will insert additional declarations immediately before the previous line. 
#endif // ! defined (AFX_LUPA_H 24CE8B94 J7D8F_11D2_958FJ)00000000000 INCLUDED_) 
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// Lupa. cpp : implementation file 

^include "stdafx.h" 
# include . n DCM . h " 
^include "Lupa.h" 

#ifdef _DEBUG 
#define new DEBUG_NEW 
#undef THIS_FILE 

static char THIS_FILE[] = FILE ; 

#endif 

///////////////////////////////////////////////////////////////////////////// 
// Lupa dialog 

Lupa: : Lupa (CWnd* pParent /*=NULL*/) 
: CDialog (Lupa : : IDD, pParent) 

{ 

//{ { AFX_DATA_INIT (Lupa ) 

l_height = 129*theApp. app_Resolutio"nScaleFactor; 

l_width = l_height; 

l_zoom =2.0; 

l_optimize=f alse ; 

//} }AFX_DATA_INIT 

l_active=f alse; 

ljpreactive^f alse; 

1_DC=NULL; 

1_DIB=NULL; 

l_Data=NULL; 

this-Mni-tialize (CSize (l_width, l_height) , 1.0, l_zoom) ; 

} 



Lupa : : -Lupa ( ) 
{ 

if(l_DC ! th«App.app_Metheus) delete 1_DC; 

if(l_Data) delete [] lJData; 

if(l_DIB) 

{ 

MetheusDeleteCompatibleBitmap ( l_DC->GetSaf eHdc ( ) , 1_DIB) ; 

} 

} 

void Lupa : : DoDataExchange (CDataExchange* p.DX) 
{ 

CDialog: : DoDataExchange (pDX) ; 
/ / { { AFX_DATA_MAP { Lupa ) 

DDX_Text (pDX> I DC_MAGN I FY_HE I GHT_E D I T , l_height) ; 
DDV_MinMaxLong(pDX, ljieight, 0, 1000); 
DDXJText (pDX, I DC_MAGN I FY_W I DTH_E D I T , l_width) ; 
DDV_MinMaxLong (pDX, l_width, 0, 1000); 
DDX_Text (pDX, I DC_MAGN I F Y_ZOOM_E D I T , l_zoom) ; 
DDV_MinMaxDouble (pDX, l_zoom, 0 . , 10 . ) ; 
DDX_Check(pDX, IDC_MAGNIFY_OPTIMIZE, l_optimize) ; 
//} }AFX_DATA_MAP 

} 



BEGIM_MESSAGE_MAP (Lupa, CDialog) 
/ / U AFX_MSG J1AP ( Lupa ) 

ON_EN_CHANGE { IDC_MAGNIFY_WIDTH_EDIT / OnChangeMagnif yWidthEdit ) 
ON_EN_CHANGE (IDC_MAGNIFY_HEIGHT_EDIT, OnChangeMagnif yHeightEdit) 
ON_EN_CHANGE ( IDC_MAGNIFY_ZOOM_EDIT, OnChangeMagnif yZoomEdit ) 
ON_EN_CHANGE ( IDCJ4AGNIFY_OPTIMIZE, OnChangeMagnif yOptimize) 
//} }AFX_MSG_MAP 
END MESSAGE MAP ( ) 



* LUPA initialization 
* 

* Input: 

* on-screen lupa size csScn, screen image zoom scrn_zoom, LUPA l_zoom l_img_zoom 
+ 

+ + + + + + * + * + + + + + + + + + + + + + + * * * » + + + # + + + 

void Lupa: : Initialize (CSize csScn, double scrn_zoom, double iupa_zoom) 
{ 

l_scnSize=CSize (csScn.cx, csScn.cy) ; 

l_zoom=lupa_zoom; 

l_img_zoom= i scrn_zoom; 

double ivzf=1.0 / (scrn_zoom*lupa_zoom) ; //combined inversed zoom 
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l_imgSize=CSize ( (long) (1+ivz. _scnSize .cx) , 

(long) (1+ivzf *• l_scnSize . cy) ) ; 
l_dragRect=CRect (0, 0, 0, 0) ; 

} 
* 

* Lupa message handlers 
* 

+ # + + + * + + * + «^ + #» + +* + + + ** + ** + + + * + + + ** + *** + + ** + + * + + + ** + +*+ + + + * + ***+ + **-*r + / 

BOOL Lupa : : OnlnitDialog () 
{ 

CDialog: : OnlnitDialog ( ) ; 
l_width=l_scnSize. cx; 
l_height=l_scnSize . cy; 
UpdateData (FALSE) ; 

return TRUE; // return TRUE unless you set the focus to a control 
// EXCEPTION: OCX Property Pages should return FALSE 

> 

void Lupa: :OnChangeMagnif yWidthEdit ( ) 

{ * 
UpdateData (TRUE) ; 
if (l_width<10) l_width=10-; 
l_scnSize . cx=l_width; 

Initialize (l_scnSize,. l_img_zoom, l_zoom) ; 

void Lupa :: OnChangeMagnif yHeightEdit t ) 

( • 
UpdateData (TRUE) ; * 
if (l_height<10) l^height^lO ; 
l_scnSize . cy=l_height ; 

Initialize ( l_scnSize, l_img_zoom, l_zoom) ; 

J 

void Lupa: : OnChangeMagnif vZoomEdit () 
{ • 

UpdateData (TRUE) ; 

if (l_zoora<l . 5) l_zoom=1.5; 

Initialize ( l_scnSize, l_img_zoom, l_zoom) ; 

\ ■ 

void Lupa : : OnChangeMagnif yOptimize < ) 
{ 

UpdateData (TRUE) ; 
l_optimize= ! l_optimize; 

} 



y + * + + + + + + * + + + + + » + + + + + + + + + + + * + + + + + + + + * 

* 

* Move lupa over the image CDC, responding to (dis) activated status 
* 

void Lupa : :Move {CPoint& a, CPoint&'rel, Image* pBmp, CDC* pDC) 
( 

CPoint p; 
CRect rO,rl; 
CSize da=a-rel; 

if(l_active) // active lupa was requested 
( 

// Remove lupa resizing rectangle*, if any 

if (l_dragRect. bottom != 0) 

{ 

Draw (pDC) ; 

l_dragRect=CRect (0, 0, 0, 0) ; 

} ' 



if ( ! i_preactive) //was not active before 
{ 

r0=CRect (0, 0, pBmp->rtr sraS . crScreen .Width ( ) , pBmp->m_smS . crScreen . Height ( ) ) ; 
Reset_l_DC(pDC, rO) ; 

r0=CRect (0, 2, 0, 2) ; // dummy update area 

) - 
else rO=l_scnRect ; // Remember previous image area 
// Compute new lupa zoom area 
rl=SetScrnRect (a) ; 

7/ For Metheus: ignore rectangles not fully inside the image area 

if ( theApp.app_Metheus ! pBmp->m_smS . Sc'reen_in_Image ( r 1 , 3 ) ) return; 

l_scnRect=rl ; 

// Compute and redraw update" rectangles 
bool bu=Update2 (rO, rl) ; 
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if {theApp. aDp_.Me theus) 
{ 

rO. InflateRect (2, 2) ; 

pBmp->DisplayDIB(pDC, rO, pBmp->m_smS . Screen_to_Image (rO) , CPoint (da) ) ; 
/* 

MetheusLoadlmageFromDIB (pDC->GetSaf eHdc () , 1_DIB, theApp. app_DynamicPaletteStart, 

rO . left, rO . top, rO . Width { ) , rO . Height { ) , 
rO.left, rO.top) ; 

*/ 

} 

else 
i 

rO.OffsetRect {-da) ; 

pDC->BitBlt ( rO . lef t , rO . top, rO .Width { ) , rO . Height ( ) , 1_DC, 
rO . lef t , rO . top, SRCCOPY) ; 

} 

if(bu) 
{ 

if (theApp. app_Metheus) 
{ 

rl. InflateRect (2, 2) ; 

pBmp->DisplayDIB (pDC, rl, pBmp->m_smS . Screen_to_Image {rl ) ,CPoint {da) ) ; 
/* 

MetheusLoadlmageFromDIB (pDC->GetSaf eHdc ( ) , 1_DIB, theApp. app_DynamicPaletteStart , 

rl . lef t, rl . top, rl .Width (), rl .Height () , 
rl . left, rl . top) ; 

*/ 

} 

else * 
{ 

rl .Of fsetRect (-da) ; 

pDC->BitBlt (rl .left, rl . top, rl .Width {) , r 1. Height () ,1_DC, 
rl . left, rl . top, SRCCOPY) ; 

} 

} 

// Zoom image into new area 

rO=l_scnRect ; 

rO .Off setRect {-da) ; 

CRect ir^Set ImgRect (pBmp->m_smS . Screen_to_Image (a) ) ; 

if (l_optimize) 

{ 

Image* kadr-new Image {) ; 

if { ! kadr->Initialize{3 + ( (7+ir. Width ( ) ) /8) ,8* ( (7 + ir .Height () ) / 3) , 

pBmp->m_Bytes_per_Pixel) ) ' 

I 

delete kadr; ' 

pBmp->DisplayDIB tpDC, rO, ir , CPoint { 0, 0) , false) ; 
pDC->DrawEdge (& (rO) , EDGE_BUMP, BF_RECT) ; 
return; 

> 

kadr->m_pPal->p_actire=false; // no palettes for lupa ! 
pBmp->GetSubimage {kadr, ir . left, ir . top) ; 
kadr->TR_HistStretch (1, false) ; 

kadr->DisplayDIB(pDC, rO, CRect {0, 0, ir .Width () -1', ir . Height {) -1 ) , CPoint {0, 0) , false) ; 
delete kadr; 

) 

else 

{ • 

pBmp~>DisplayDIB(pDC, rO, ir, CPoint (0, 0) , false) ; 

) 

pDC->DrawEdge (& (rO) , EDGE_BUMP, BF_REGT) ; 

} 

else // disactivated lupa was requested 

{ - ■ 

if (l_preactive) ~" 

p=l_scnRect . TopLeftO -da; 
if (theApp. app_Me theus) 

pBmp->DisplayDIB (pDC, l_scnRect , pBmp->m_smS . Screen_to_Image {l_scnRect) ,da) ; 
/* 

MetheusLoadlmageFromDIB (pDC->GetSaf eHdc () , 1_DIB, theApp, app_DynamicPaletteStart , 

p.x, p.y, l_scnRect .Width () , l_scnRect. Height ( ) , 
p.x,p.y);. 

*/ 

} 

else . 
{ 

■ pDC->BitBlt (p.x, p. y, l_scnRect . Width ( ) , l_scnRect . Height ( ) , 1_DC, 
p. x, p. y, SRCCOPY) ; 

delete 1 DC; 
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} 

l_preact ive«f alse ; 
1_DC=0; 

} 

} 



* 

* Represents update region as two rectangles 
+ Returns false if only one rectangle "a" must be updated 
or true if both "a" and "b" 

+ 

bool Lupa: :Update2 (CRect & a, CRect & b) 
{ 

if (a.WidthO !=b.Width() II a . Height ()! =»b . Height ( ) ) return false; 

if (a.EqualRect (&b) ) // coinciding rectangles 

{ 

a=CRect (a . left, a. top, a . left, a. top) ; 
return false; 

} 

CRect al; al .CopyRect ( Sa) ; 
CRect bl; bl .CopyRect (&b) ; 
if (a.leftob.left && b.leftoa. right) 

< 

if (a. top<=b. top b. top<=a. bottom) 
{ 

al-CRect ia . left, a . top, b. left , a .bottom) ; 
bl=CRect (b. left, a. top, a. right, b. top) ; 

) 

else if (a . top<=b. bottom b. bottom<=a .bottom) 
( 

al=CRect (a. left, a . top, b. left, a .bottom) ; 
bl=CRect (b. left, b. bottom, a . right, a . bottom) ; 

) 

else return false; 

} 

else if (a'. left<=b. right b . right Oa . right ) 
{ 

if (a. top<^b. top b. top<=a. bottom) 
{ 

al=CRect (a. left, a. top, b. right, b. top) ; 
bl=CRect (b. right, a . top, a . right, a .bottom) ; 

) 

else if {a . topob .bottom && b . bottom<=a . bottom) 
{ 

al=CRect ( a . lef t , b . bottom, b . right , a . bottom) ; 
bl=CRect <b. right, a. top, a. right, a. bottom) ; 

} 

else return false; 

} 

else return false; 

a . CopyRect (&al) ; 

b. CopyRect (&bl) ; 
return true; 

} 



* 

* Copies current screen image-- into* lupa CDC 1_DC 

void Lupa: :Reset_l_DC (CDC * pDC, CRectfi scrolled_client) 
{ 

int w=scrolled_client .Width {) ; 

int h=scrolled_client . Height ( ) ; 

l_preactive= s true; 

if ( theApp.app_Metheus) 

< 

l_DC=pDC; return; 

w=2048; h-2560; 
// BYTE buffer 

if(l_Data) { delete [] l_Data; l_Data=NULL; } 

if (1_DIB) 

{ 

MetheusDeleteCompatibleBitmap (pDC->GetSaf eHdc < ) , 1_DIB) ; 
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1 DIB=NULL; 

} 

l_DIB=MetheusCreateCompatibleBitmap (pDC->GetSafeHdc {) , w,h) ; 

l_Data=new BYTE[w*h*3]; 
HOC hdc=pDC->GetSaf eHdc ( ) ; 

BOOL bl=MetheusGetImageIntoData (hdc, 1_DIB, theApp.app_DynaiaicPaletteStart, 

(UCHAR*) l_Data, w, h, 2*w, 16, 
0, 0, w,h, 
0,0) ; 

if (bl — FALSE) 
{ 

Beep (700, 150) ; 

AfxMessageBox ("MetheusGetlmagelntoData Failed") ; 

} 

BOOL b2=MetheusLoadImageFromData (hdc, 1_DIB, theApp . app_DynamicPa let teS tart , 

l_Data, w, h, 2*w, 16, 
0, 0, w,h, 

scrolled_client . lef t, scrolled_client . top) ; 

if (b2==FALSE) 
{ 

Beep{700, 250) ; 

AfxMessageBox { "MetheusLoadlmageFrcmData Failed") ; 

} 

} 

else 
{ 

if(l_DC) { 1 DC->DeleteDC () ; delete 1_DC; 1_DC=0; } 
• i_DC=new CDC (7; 
l_DC->CreateCompatibleDC (pDC) ; 
CBitmap b; 

b.CreateCompatibleBicmap(pDC, w, h) ; 
l_DC->SelectObject (&b) ; 

l_DC->BitBlt (0, 0, w, h, pDC, scrolled_client . left, scrolled_client . top, SRCCOPY) ; 
b. DeleteObject () ; 

} 

} 
* 

+ Interactively resize the Lupa region on the image 

** + *** + + + *.** + ** + + + + + + + * + * + + + *** + + + + ** + * + *^ 

bool Lupa : : Resize (CDC *pDC, CPoint center, CPoint vertex) 

{ 

Draw(pDC); // remove old rectangle 
CPoint z=vertex-center; 
int a~<abs(z.x) )<<1; if(a<32) a~32; 
int b=(abs(z.y) )«1; if(b<32) b=32; 
Initialize (CSize (a, b) , l_img_zoom, l_zoom) ; 
l_dragRect=SetScrnRect (center) ; 
Draw(pDC); // draw new rectangle 
return true; 

} 

bool Lupa : : Resize (CDC *pDC> 
{ 

Draw(pDC); // just remove old rectangle 
l_dragRect=CRect (0,0,0,0); 
return true; 

} 
* 

* Output Lupa parameters" into a. string (used in status bar) 

+#++**++***++*+++*+*+**+**+^***+******+*+*+* + ** 

CString Lupa : : toString ( ) 

{ 

CString s; 

s. Format ("Screen size sdx^d, zoom=% . 21f ", l_scnSize . cx, l_scnSize. cy, l_zoom) ; 
return s; 

} 
* 

* Draw Lupa region rectangle 

void Lupa: : Draw (CDC *pDC) 
{ 

int dmode=SetROP2(pDC->m_hDC, R2_NOT) ; 
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pDC->MoveTo ( l_dragRect . TopLe . ) ) ; 
pDC->LineTo ( l_dragRect . right , l_dragRect . top) ; 
pDC->LineTo ( l^dragRect . right , l_dragRect . bottom) 
pDC->LineTo (l_dragRect . left, l_dragRect . bottom) ; 
pDC->LineTo(l_dragRect .TopLeftO ) ; 
SetR0P2 <pDC->m_hDC, dmode) ; 



